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ABSTRACT 

We estimated the dynamical surface mass density £ at the solar position between 
Z=1.5 and 4 kpc from the Galactic plane, as inferred from the kinematics of thick disk 
stars. The formulation is exact within the limit of validity of a few basic assumptions. 
The resulting trend of matches the expectations of visible mass alone, and no dark 
component is required to account for the observations. We extrapolate a dark matter 
(DM) density in the solar neighborhood of Oil 111M0 pc~ 3 , and all the current models 
of a spherical DM halo are excluded at a confidence level higher than 4<t. A detailed 
analysis reveals that a small amount of DM is allowed in the volume under study by the 
change of some input parameter or hypothesis, but not enough to match the expecta- 
tions of the models, except under an exotic combination of non-standard assumptions. 
Identical results are obtained when repeating the calculation with kinematical measure- 
ments available in the literature. We demonstrate that a DM halo would be detected by 
our method, and therefore the results have no straightforward interpretation. Only the 
presence of a highly prolate (flattening q >2) DM halo can be reconciled with the obser- 
vations, but this is highly unlikely in ACDM models. The results challenge the current 
understanding of the spatial distribution and nature of the Galactic DM. In particular, 
our results may indicate that any direct DM detection experiment is doomed to fail, if 
the local density of the target particles is negligible. 
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1. INTRODUCTION 



Today, it is widely accepted that a dominant fraction of the mass in the universe is in the form of 
a non-luminous (dark) matter (DM), whose nature is still unknown. Despite this general agreement, 
and decades of investigations, its spatial distribution and its density in the solar neighborhood 
(Pq,dm) are still poorly constrained by the observations. The local densitj|§ extrapolated from 
the Milky Way rotation curve and other observationa l data, assuming a spherical DM halo, is 
in the range 5-13 mM pc -3 (IWeber fc de Boerl I201C1 ). However, /9©,dm could be significantly 
different in the case of an oblate or prolate halo, or in the presence of a DM disk. Our poor 
knowledge of the DM density distribution is very unfortunate, since this information is crucial 
to clarify the nature and properties of DM itself. The shape of the dark halo, quantified by 
its shortest-to-longest axis ratio q, is in fact an important diagnost ic on the natu re of the DM: 
for example, a round halo (q ~1) is expected by hot DM models ([Peebles! Il993l ). while a very 
fiat halo (q ~0.2) is pr eferred if cold molecu l ar gas or ma ssive decaying neutrinos are the main 
constituent of the DM (jPfenniger et al.l 11994 : ISciamal Il990l ) . Noticeably, great effort is currently 
spent in experiments for direct DM detection. The results of these experiments are degenerate 
betw een the unknown interaction cros s-section of the searched particles and their local density 



(e.g., iGaitskelll 12004 : lAprile et al.ll2005l ). Thus, most of the works estimating the properties of the 



Weakly Interacting Massive Particles (WIMPs) hav e so far drawn their c onclusions assuming the 
local DM density of the Standard Halo Model (SHM. I Jungman et al.lll996l ). Pq,dm = 8 mM pc~ 3 . 
Clearly, observational constraints on the DM density at the solar position and on the flattening of 
the dark halo are key components for revealing the secrets of the Galactic DM. 

"Weighing" the Galactic disk by means of the spatial distri bution and kinematics of its stellar 
component is an ancient art, dating back nearly one century (IKapteynl Il922l : lOorti Il932l ). The 
difference between the measured mass and the visible mass provides an estimate of the amount of 
DM in the volume under analysis, and constraints on the shape of the DM halo can be derived. 
The fundamental basis for this classical measurement is the application of the Poisson-Boltzmann 
and Jeans equations to a virialized system in steady state. This allows us to estimate either the 
local density at the solar position or the surface density (mass per unit area) of the mass within 



1 Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemi- 
sphere, Chile (proposal IDs 075.B-0459(A),077.B-0348(A)) 

2 This paper includes data gathered with the 6.5-meter Magellan and the duPont Telescopes, located at Las 
Campanas Observatory, Chile. 

3 The DM mass will be given in mM throughout the paper, where 1 mM = 10~ 3 M© = 0.038 GeV. 
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a given volu me. iGarrido Pestana & Eck hardtl (|2010i ) recently argued against the reliability of this 
method, but ISanchez-Salcedo et al.l (|201ll ) demonstrated that it can be applied to the Galactic 
disk, because it is in equilibrium with the Galactic potential. The measurements of the dynamical 
mass in the solar ne ighborhood are abundant in the literature, and all but few works (|Bahcarllll984l ; 
Bahcall et al.lll992l ) came to the same overall conclusion that "there is no evidence for a significant 
amount of disk DM 1 . Different interpretations of this statement have been presented though: while 
the measurements of the surfa ce density usually match the expectations of the visible mass plus 
a classical spherical DM halo (IKuiiken fc Gilmore! ll989l; Flvnn Fuchj Il994j : Isiebert et alJhood ; 
Holmberg &: Flvnn 2004 ; Bienavme et al. 20061). the estimates of the local volume density usually 



find a much lower quantity of DM (IKu iiken fc Gilmore 



2000 



Korchagin et al. 



2003 



1989 



Creze et al 



1998 



Holmberg &: Flvnn 



de Jong et al.ll2010i ) , although they are compatible with the presence of 



a classical halo within the uncertainties. Despite many decades of investigation, little progress has 
been made beyond this statement: for example, the most recent measurements are still compatible 
with both the expectation s of the SHM and the complete absence of DM in the solar neighborhood 



Garbari et al 



2011 



). Some weak constra ints on the properties of the DM halo were given by 



Creze et al. (jl998l ) and Bienavme et al. ( 2006 ). who claimed that a DM halo flattening q < 0.5 is 
excluded by the observations. 

The main limitations to the measurements of the dynamical mass are imposed by the great 
observational effort required, because the information about the kinematics and spatial distribution 
of a vast number of stars is needed. Some approximations were thus always introduced in the 
formulation, accurate on ly at low Galactic h eights , whose valid i ty ha s since been questioned in the 
literature. For example. ISiebert et al.l (120081 ) and ISmith et al.l (|2009l ) found that at only Z=l kpc 
from the Galactic plane the potential is not separable in the radial and vertical co ordinates, and 



neglec ting the non-diagonal term of the dispersion matrix could introduce a bias. iGarbari et al 



(|201ll ) found that this beco mes an issue at onl y Z=0 . 5 kpc , and argued against the use of the 
approximate formulation of iHolmberg &; Flvnnl (|2000l . |2004) at high Z. As a consequence, all 
previous investigations were limited to Z < 1.1 kpc from the Galactic plane. The amount of DM 
expected in this volume is smal l compared to ob s ervat ional uncertainties, and no strong conclusion 
could be drawn. Moreover, onlv lKorchagin et al.l (|2003l ) directly calculated the mass density from an 
analytical expression, while the other investigations estimated it by comparison of the observational 
quantities with the expectations of a Galactic mass model. 

In this paper, we propose a new approach to estimate the surface mass density of the Galactic 
disk. Our formulation leads to an analytical expression that is exact within the limits of validity 
of a few basic assumptions. The surface density can thus be directly calculated at any distance 
from the Galactic plane. The calculation requires a knowledge of the spatial variations of the 
three-dimensional kinematics of a test stellar population. Our calc ulations are based on the results 



of ou r recent investigation of the Galacti c thick disk kinematics (jCarraro et al 



2005 



Moni Bidin 



2009i ). Similar results were presented by lMoni Bidin et al.l (|2010l ). who searched for a signature of 
the dark matter disk, predicted by the merging scenario of thick disk formation. 
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2. THE THEORY 

In the following, we will use the cylindrical Galactic coordinates (R,9,Z), where R is the 
Galactocentric distance, 9 is directed in the sense of Galactic rotation, and Z is positive toward 
the north Galactic pole. The respective velocity components are (R, 9, Z)=(U, V, W). 

Our formulation is based on the integrated Poisson equation in cylindrical coordinates 

- 4nGZ(Z) = ±JL(RF R )dz + 2 • [F Z (Z) - F z (0)i (1) 

where G is the gravitational constant, is the surface density of the mass comprised between 

±Z from the Galactic plane, and the radial and vertical component of the force per unit mass, Fr 
and Fz respectively, can be expressed through the Jeans equations 



d<j> = 1 d( P U*) 1 d(pUW) U2- V 2 i d{pU) 
R dR p dR p dZ R p dt 1 U 



F z = 

z dz 



d(j) _ 1 \d{pW 2 ) pUW d(pUW) d(pW) 



(3) 



dZ R dR dt 

where p is the volume density, and is the gravitational potential. The following basic assumptions 
can be made, for symmetry reasons, for the stellar population under study: 

I The test population is a virialized system in steady state. All the temporal derivatives are 

therefore null: 

d{ P U) _ d{ P W) 



dt dt 



0. (4) 



II The vertical component of the force per unit mass is null on the plane: 

F z (0) = 0. (5) 

III The trend of velocity dispersions with Z is symmetric with respect to the plane: 

Ui{Z) = <n(-Z), (6) 

with i = U,V, W. 



IV UW is antisymmetric with respect to the Galactic plane: 



UW(Z) = -UW(-Z). (7) 
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V The net radial and vertical bulk motion of the trace population is null: 

U = W = 0. 



(8) 



Sanchez-Salcedo et al.l (120111) demonstrated the validity of a ssumption (JT|), that was recently ques- 
tioned in the literature (jGarrido Pestaha fc EckhardtJl2010l ). The hypothesis (|TTJ) , (IIIip . and (|IVp 
are required for symmetry reasons, and they imply that the integral of the velocity dispersions in 
±Z is twice the integral between zero and Z, while the integral of UW is null, and UW(0)=0. It 
can be easily seen that assumption (|IV[) directly results from the Jeans e quations, if the symmetry 
requirements Fr(Z)=Fr(— Z) and Fz(Z) = —Fz(—Z) are to be satisfied. iMoni Bidin et al.l (|2010l ) 
also showed that the cross-term must be assumed anti-symmetric with respect to Z, else the cal- 
culation leads to unphysical results. Recent investigations de tected a non-null mean radial motion 



of stars ( Casetti-Dinescu et al 



2011 



Moni Bidin et al.ll2012l ). However, this is not a consequence 



of a global motion of disk stars, but it is rather a signature of local kinematical substructures, so 
that assumption (JVJ) is still valid as a general property of disk stars. 



Inserting Equations ^) and ([3D into Equation (pTJ) , and making use of the assumptions (J]) 
to (jVJ), the resulting equation can be analytically solved for at any R and Z, if the three- 

dimensional analytical expressions for the kinematics and the mass density distribution of the test 
population are given. To provide this extensive input, we make use of observational results only. 
The spatial distribution of the test population is fixed by the following assumptions: 

VI The volume density of the test population decays exponentially both in the radial and vertical 
direction, with exponential scale height and length hz.p and iiR, p , respectively: 



p(R, z) = p • exp 



Z — Zq R — Ra 







hz, P hn tP 
where p & is the density at (R, Z) = (Rq, Zq). 
VII The scale length of the test population, h/j )/9 , is invariant with respect to Galactic height: 

dh RjP 



0) 



dZ 



0. 



(10) 



These hypotheses are a classical re presentation of Galact i c disk -like populations, and (|VIip is con- 
firmed by the empirical results of ICabrera-Lavers et al.l (|2007l ) . Theoretical considerations sug- 
gest tha t the vertica l dens ity profile should be much closer to a sech 2 function than an expo- 
nential (jCamml Il950l . [1952J) , but observational evidences confirm that the Galactic thick disk, 
object of our investigation, is well described by the double exponential law of Equation ([9]) (e.g., 
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Hammerslev et al.lll999l ; ISiegel et al.l 120021 ; I Juric et al.ll2008l ). The assumed vertical decay is surely 
an accurate fit of the disk density at our large Galactic heights, where even the sech 2 profile ap- 
proximates to an exponential decline. 

The calculation is simplified by two additional working hypothesis: 



VIII The rotation curve is locally flat in the volume under study: 

dV 



(ii) 



IX The disk test population is not flared in the volume of interest: 

dh z , p 



dR 



0. 



(12) 



These are very common assumptions, whose function in our approach is neglecting few small terms, 
which complicate the formulation while introducing only second-order corrections. The observa- 
tional evidences in their support will be discussed in Section 14.21 and I4.4[ where the influence of 
their break-down will also be analyzed. 

It is easily demonstrated that, by means of our assumptions, the final expression for the surface 
density reduces to 



2vrGS(i?,Z) 



dZ 



a 



hz, P 
1 

~R 



a 



V 



■dz + 



o RhR, p 
z da 2 v , / 2 



-dW dz + 



2 
R 



1 



l 



+2 



dUW 
dR 



h R)P J dR 
h z p Jo oR 



This equation is exact within the limits of validity of the underlying assumptions, and it can be used 
to calculate the surface density of the mass within ±Z from the Galactic plane at the Galactocentric 
distance R. The estimate requires the knowledge of the two parameters hz.p and h^, the vertical 
trend of the three dispersions and of UW, and the radial derivative of a\j(Z), ay(Z), and UW(Z). 
A great observational effort, impossible until the last decade, is required to collect this quantity of 
information. Nevertheless, this will be eventually gat hered by modern extensive sur veys such as 
the Sloan Digital Sky Survey (SDSS; 



Synoptic Survey Telescope (LSST; llvezic et al 



York et al 



2000), GAIA ([Wilkinson et al J 120051 ). and Large 
2008). The measurement of the mass distribution 



in the Galaxy will be then possible with unprecedented detail. 

As shown in Section [3j the data used in the present study provide no information about the 
radial behavior of the kinematics, and we are forced to model it with an additional assumption. 
We adopt the most likely hypothesis: 
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X UW and the square of the dispersions exponentially decay in the radial direction, with the 
same scale length as that of the volume mass density. 



The effects of alternative radial trends will be analyzed in Section 14. li Assumption (|X|) is observa- 



tion H3J 
1982h . ar 



tionally confirmed only for cr\y (jvan der Kruit Sz Searldll98ll . I1982 ). and its extension to the other 
components relies on the assumption of radially constant anisotropy, i.e. the radial constancy 
of the ratio of t he dispersions. However, bot h theoretical calculations and observations support 
this hypothesis: ICuddeford &: Amendtl (|1992l ) demonstrated, by means of numerical integration 
of orbits, that it is the best representation of the radial trend of dispersions for R < 9 kpc, and 
the most deta i led ob servations available to date also co nfirm it (ILewis Freemanlll989i ). although 
Neese Yossl (| 19881 ) prefer a linear radial decay for a\j. iMoni Bidin et al.l (|2012l ) also find that the 
observational data are consistent with assumption (|X)) . 



Using assumption (jX}, Equation [13] can be simplified substantially, and after a calculation 
involving only simple integrals and derivatives, we obtain 



E(Z) 



1 



2ttG 



a 



h- I a(jdz + k 2 - I aydz + k 3 ■ UW + 

JO h Z,p 



da 



w 



dZ 



where 



k 2 
k 3 



R Q ■ h R:P h? R p ' 
1 



3 2 



(14) 

(15) 
(16) 
(17) 



3. RESULTS 



Our estimate of the surface mass density is based on the results of IMoni Bidin et al.l (|2012l ). 
who measured the kinematics of the Galactic thick disk, and its variation with distance from the 
plane, between Z=1.5 and 4.5 kpc. In brief, they studied a sample of ~400 thick disk red giants 
toward the south Galactic pole, vertically distributed with respect to t he Galactic plane, and derived 
their three-dimensional ki nematics by means of 2MASS photometry fSkrutsk ie et al.ll2006l ). SPM3 
absolute proper motions (jGirard et al.ll2004l ). and radial velocities ()Moni Bidinl 120091 ) . Thus, the 
data do not provide any information about the variation of the kinematics with Galactocentric 
distance, that will be modeled with a ssumption dXl). The sur face mass density will therefore be 
calculated by means of Equation (|14p . IMoni Bidin et al.l (|2012l ) detected a clear increment with Z 
of all the dispersions, well represented by a linear fit. We wi l l ado pt for a\j(Z), ay(Z), and ay^(Z) 



the relations given by equations (3)-(5) of IMoni Bidin et al.l (120121 ) . The vertical trend of UW will 
be taken from the linear fit shown in Figure 2 of lMoni Bidin et al.l (|2010l ) , which yields the solution 

TJW = (1522 ± 100) + (366 ± 30) • [Z - 2.5) 



km 2 s 2 , 



(18) 
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2 3 4 

Z [kpc] 

Fig. 1. — Observational results for the absolute (upper panel) and incremental (lower panel) surface 
mass density, as a function of distance from the Galactic plane (black curves), compared to the 
expectations of the models discussed in the text (thick grey curves). The dotted and dashed lines 
indicate the observational la and 3<r strip, respectively. 
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where Z is in kpc. 

Measuring the kinematical properties and the spatial distribution of the test population by 
means of the same obse rved stars would b e highly desirable. Unfortunately, our sample is not 
suitable for this: even if Girard et al. (j20od ) used it to estimate the thick disk scale height, their 
measu rement could have been affected by the thin disk contamination unaccounted for in their 
study (jMoni Bidin et al,ll2012l ). However, our sample can be considered representative of the thick 
disk kinematics at the solar position, whose spatial distribution has been extensively studied in 
the last three decades. This would not be true if the thick disk was a mixture of sub-populations 
with different properties. In this case, the spatial distribution of the sub-population under study, 
or the dominant one, should be preferred. Nevetheless, to our knowledge, there is no evidence in 
literature supporting such high degree of inhomogeneity in the thick disk, nor that the sample under 
study is for some reason peculiar. Finally, we fix R = R & = 8.0±0.3 kpc, and h.R iP = 3.8±0.2 kpc, 
h^ p = 900±80 pc from the average of sixteen and twenty-one literature measurements respectiveljfl. 
The impact of these parameters o n the final results will be analyzed in Section 14.11 The errors were 
defined, as in iMoni Bidin et al.l (|2010l ). from the statistical error-on-the-mean. If the literature 
estimates can be considered independent measurements of an underlying quantity, this must be the 
most rigorous estimate of the true uncertainty. However, the scale height and length of the Galactic 
thick disk are traditionally poorl y constrained, although t he measurements converged considerably 
in the last years. For this reason. IMoni Bidin et al.l (|2010l ) also considered enhanced errors (0.4 and 
0.12 kpc for 1ir iP and h^ jP , respectively), possibly more representative of the true uncertainties. In 
this case, the final errors on are enhanced by about 50%, and the significance of the results 

thus decreases by a factor of ~0.7. Nevertheless, adopting these larger errors would not affect the 
general conclusions of our work noticeably because, as will appear clearer later, the significance of 
the results remains very high even in this case. 

The error on S(Z) was calculated from the propagation of the uncertainties on the four kine- 
matical quantities and the three parameters. The resulting equation is particularly cumbersome, 
but it can be obtained through simple derivations of the terms on the right hand side of Equa- 
tion (I14j) . Monte-Carlo simulations, kindly provided us by the referee, confirmed that the final 
error is a good estimate of the uncertainty propagated from the ones associated to the quantities 
entering in the calculations. The simulations were performed repeating the estimate of after 
varying the input quantities, each one randomly drawn from a Gaussian distribution with mean 
and dispersion given by the assumed value and error, respectively. 

The resulting vertical profile of the surface mass density ~S{Z), is shown in the upper panel 



featnatunga , et al.l (Il989l); lYamagata fc Yoshiil Jl992l); Ivon Hippel fc Bothunl (I l993); 



(199. 



5) 



Larsenf (Il996l); iRobin et al.l (Il996l); ISpagna et al 



Chen et al 



Oiha et al.l dl999h: Ichiba fc Beers! J200oll 
Larse^T^TliumpIirevsl j2003l);ICabrera-Lavers et al 



iOOahlGirard et al 



( 1996) ; iNg et all ( 1997) 



2008) ; ICabrera-Lavers et al.l (120071'); lArnatoddir et al 



Beers fc Sommer-Larsenl 
Buser et al.l Jl99£ 1: 



(hoOlh; iLopez-Corredoira et al.l (12003); Isiegel et al.l J2OO2I) 



(2 



(2008) 



20081 ); Ide Jong et all |2010l ); ljust et alj l|201lh ; Idiang et al.1 (|201ll ) 



200d);IVaHenari et al.1 12009)'; I Brown et al 
Bilii' (it al.1 J2OO8I ); ljuric et all j|200sl ): IVeltz et al 
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of Figure [TJ We find £(1.5 kpc)=55.6±4.7 M pc~ 2 , and the profile is nearly flat, increasing by 
only ~2 Mq pc~ 2 up to Z=3 kpc and eventually bending down, while the error monotonically 
increases up to 10 M Q pc~ 2 at Z=A kpc. The decrease of T*(Z) in the last kpc is unphysical, but it 
is negligible compared to the errors, corresponding to only 0.07o\ This could be corrected by a tiny 
change of the input parameters, for example increasing h/j p by 0.1 kpc, but we will not manipulate 
them to avoid the introduction of arbitrariness in the results. If the solution is extrapolated to 
Z=0, the unphysical result £(0) ^ is obtained. This is due to the fact that the vertical trend 
of the kinematical quantities was derived from a linear fit. This is a good approximation in the 
.Z-range under analysis, but the dispersions should depart from it at lower heights, bending down 
with steeper gradient. The most clear case is UW, which should be zero on the plane for symmetry 
reasons, but the linear fit returns UW(0) = 607 km 2 s -2 . The extrapolation of the results thus 
fails to account for the changed gradient of the kinematical quantities and, as a consequence, of 
£(Z). In fact, as shown later (Figure [6j), a steeper curve of the surface mass density at lower Z 
(and a lower extrapolation t o Z=0) is recovered when using the steeper gradients of kinematics of 



Casetti-Dinescu et al.l (|201ll ). It can be noted that the linear fit is a worse approximation in this 
case, as reflected by T>(Z) being too flat compared to any model, because the dispersions vary their 
gradient more quickly nearer to the plane. In conclusion, the extrapolation of our results outside 
the Z-range where they were obtained is not allowed. 

In Figure [H the results are compared with the expectations for the known visible mass, in- 
dicated by the thick grey curve labeled as VIS. This was modeled as the sum of a thin layer of 
interstellar medium (ISM) of 13 M Q pc~ 2 (jHolmberg Flynnl 12000 ) . plus three stellar compo- 
nents, the halo, t he thin disk, and the thick disk, whose geometrical parameters were taken from 
Juric et al.1 (|2008l ). Their local density at Z=0 was normalized so that the total surface densit y 



of the stellar disk is £di s k(l-l kpc)=40 M pc 2 (jHolmberg Flynnll2004: iBienavme et al 



This is currently the best estimate assumed in Ga lactic mass models (e.g.. iDehnen &: Binney 



2006) 



1996 



Oiling & Merrifieldll200ll ; IWeber & de Boerll20ld ). but the local mass density of both the ISM and 
the stellar component are still affected by observational uncertainties. However, as discussed below, 
the exact value of these parameters does not have a significant influence on our results. 

The estimate of the surface mass density matches the expectation of visible mass alone, and the 
degree of overlap between the two curves is striking. There is no need for any dark component to ac- 
count for the results: the measured £(Z) implies a local DM density Pq,dm = 0±1 mM pc~ 3 . This 
estimate negligibly change s (p^.nM = 0.4zfcl.2 mM(T ) pc~ 3 ) if we assume £riisit (1-1 kpc)=35 M Q pc~ 2 



in the visible mass model (jHolmberg h Flvnnll2000l : iGarbari et al.ll201ll ). Fitting a DM halo model 
to the observations is pointless because, disregarding its exact shape, the procedure necessarily 
converges to a zero-density solution. The more complex task of building a model able to repro- 
duce our results and other observational constraints (e.g. the Galactic rotation curve, the gas disk 
flare) is beyond the scope of this paper. Nevertheless, we can compare the measurements with the 
expectations of the most popular DM halo models. 



A great quantity of models have been proposed in the literature to describe the spatial dis- 
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tribution of the Galactic DM, but only the mass between and 4 kpc from the plane at R = R@ 
is involved in the comparison with the observations. Hence, the local density Pq,dm and the halo 
flattening q are the only relevant parameters, while the exact functional form is not critical. In- 
deed, it is easy to see that different models with the same Pq.dm and q are indistinguishable in 
terms of the expected quantity of DM in the volume under analysis. We will therefore compare 
the observations with a set of selected models, chosen on the basis of the expected Pq ; dm- We will 
assume spherical models here (q=l), while the effects of varying the flattening parameter will be 
analyzed in more detail in Section 14.71 



Oiling fc Merrifieldl (|200ll ) presented a family of self-consistent models, based on a classical 
non-singular isothermal spheroid: 



Pbm(R, Z) — p c 



R,, 



Rl + Rl + (z/ q y 



(19) 



(e.g., Ivan Albada Sancisil Il986l : iKentl 1 1 98 71 ) . where the core radius R c an d the central density 



p c depend on q, so that the resulting rotation curve is independent of it (Oiling 



1995). As a 



conservative choice, we will assume, among the solutions proposed by lOlling Sz Merrifieldl (1200 ll ) 
with i?0=7.8-8.5 kpc and S(ji s k(l-1 kpc)>30 M pc -2 , the model with the minimum local density 
(Pq,dm ~10 mM pc -3 , hereafter model OM), which have i? c =8.01 kpc and p c =20.6 mM pc~ 3 . 



A more general expression of the DM halo shape is 



PDM(R, Z) = Pq,DM 



R, 



■0 



1 + 



\ Rc > 



(20) 



where the indices (a, {3, 7) characterize the radial fall-off of the density distrib ution. Various sets of 



indices have been proposed in the literature, suggesting eit her "cuspy" (e.g.. [Navarro et al 



Moore et al. 


1999 


Binnev & EvanskoOl: 


Ludlow et al.ll 2009) 


2005; 


Naravan et al. 


2005: 


Gentile et al. 


2007; 


Salucci et al. 



resentative models were analyzed by 



1997; 



de Boer et al 



Weber fc de Boer 



who fixed the best-fit parameters 



on the basis of the most recent observational constraints. Among these, we will c onsid er the model 
which better fits the observations (in terms of the lowest x 2 )> the lNavarro et al.l ( 19971 . NFW) pro- 
file (a = (3 = 1, 7 = 2) with i? c =10.8 kpc and Pq : dm=8A mM Q pc -3 . This value coincides with the 



Stan dard Halo Model density usually assumed in direct DM detection experiments (jJungman et al 



1996), and it will be referred to as SHM. We will also consider the NFW profile with the lowest 
local density (R c =20 kpc, /0©,DAf=6.1 mM pc~ 3 , hereafter model N97), and the model with the 



minimum local DM density, a pseudo-isothermal profile (a = 0, /3 = 2, 7 = 1, Ide Boer et al.ll2005l . 



hereafter model MIN) with R c =5 kpc an d /3^ i nM=5.3 mMp pc 3 . This value is usua lly assumed 



as a lower limit for the local DM density (IGarbari et al 



2011 



Weber fc de Boerll2010h . 



In Figure [H the measured vertical profile of the surface density is compared to the expectations 
of the selected DM halo models. The OM model is excluded at the 8<r level, while the SHM at the 
6<7 level. Even the MIN model, with the minimum local density extrapolated from the Galactic 
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rotation curve (jGarbari et al.ll201ll ). is A.la more massive than the observed curve. If we assume 
^disk(l-l kpc)=35 M pc~ 2 in the visible mass model, the results have a lower but still very 
significant departure: even the MIN model would still be excluded at the 3. 6a level. 

It could be argued that integrating the linear vertical profiles of the velocity dispersions in the 
range Z=0-4 kpc implicitly assumes their extrapolation down to Z=0, while they were measured 
only for Z >1.5 kpc. Indeed, the vertical trend of the dispersions is not expected to be strictly 
linear, but to be shallower at higher distance from the plane (see, for example, the discussion of 



Caset ti-Dine scu et al.ll2011). Ne vertheless, assuming for Z <1.5 kpc the steeper relations measured 
by ICasetti-Dinescu et al.l (|201ll ). the resulting surface density decreases by only 0.5 M© pc~ 2 . 



Arbitrarily increasing all the gradients by a factor of two for Z=0-1.5 kpc, decreases by 

5.5 Mq pc -2 , enhancing the disagreement with the DM halo models. In conclusion, the integration 
of the profiles in the entire range is likely a good approximation, and it cannot be the cause of the 
mismatch between the measurements and the model expectations. 

The increment of the surface density between 1.5 and 4 kpc (£>i.5(Z)) can be measured inte- 
grating Equations (|13p and (|14p in this interval, instead of 0-4 kpc. This provides a measurement of 
the surface density of the mass enclosed between 1.5 and 4 kpc from the plane only. This approach 
has the advantage of being free of the uncertainties related to the total amount of visible mass and 
the extrapolation of the kinematics to Z <1.5 kpc, although we have already demonstrated that 
these points do not invalidate our conclusions. In fact, all the ISM and the majority of the stellar 
mass are found below the range of integration, and a decrease of Edi S k(l-l kpc) by 5 M Q pc~ 2 
decreases the expectation of S>i.5(Z) by only 0.15 M Q pc~ 2 . 

The derived profile of T, > i^(Z) is shown in the lower panel of Figure [TJ The same conclusions 
as before can be drawn: the results perfectly match the expectation of the visible mass alone, and 
the estimated local DM density is p©,dm = 0.0 ± 0.7 mM pc -3 . Moreover, the discrepancy with 
the models comprising a DM halo is even more striking, since the OM model is 9a higher than 
the derived solution, and the SHM and MIN models are excluded at the ~ 7.5a and ~ 5a level, 
respectively. 



4. ANALYSIS 

The calculation relies on a set of kinematical measurements, ten hypotheses, and three pa- 
rameters. In this section, we will assess in more detail the reliability and robustness of the results, 
investigating if any of these input quantities can be the cause of the mismatch between observations 
and model expectations, and if the DM halo could pass undetected with our methods. Finally, al- 
ternative DM spatial distributions (non-sperical halo, dark disk, dark ring) will also be considered 
and compared to the observations. 
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4.1. The incidence of the parameters 



A different definition of the solar Galactocentric distance has a negligible impact on the results. 
Both Tj(Z) and S>j.5(4 kpc) decrease with R & , but they change by less than ±3.5 M pc -2 and 
±2.6 M pc -2 (~ 0.5cj at 4 kpc), respectively, when Rq is varied in the range 7.5-8.5 kpc. The 
precise value adopted for Rq is therefore irrelevant. 

The effects of varying the thick disk scale height and length in the range 0.6-1.2 kpc and 
3.0-4.6 kpc, respectively, are shown in Figure [2j Decreasing hz tP shifts the curve of to 
higher values, but does not affect its slope noticeably. For example, X (1.5 kpc) agrees well with 
the expectation of the most massive DM halo model (OM) when assuming hz tP =0-6 kpc, but at 
Z=A kpc they differ by ~ 4cr, because the gradient of the observed curve is too shallow. In fact, 
even assuming this extremely low scale height, the less massive model (MIN) is still 2.5cr higher 
than the curve of T l> i^(Z). On the contrary, varying the thick disk scale length affects the slope 
of T>(Z), but does not add a lot of mass to the derived solution. S>! 5 (Z) increases with hji tP , but 
even assuming h^ p =5 kpc, the MIN model is still 2a higher than both S(Z) and T l> i^(Z). 

While varying hz, P or hn iP alone is not enough to reconcile the measurements with the pres- 
ence of a classical DM halo, they can be changed simultaneously. Nevertheless, observational 
constraints prevent from obtaining a solution overlapping the expectations of DM halo models, 
because even the less massive MIN model can be roughly matched by the observations only by as- 
suming ^.^=0.65 kpc and hR iP =A.7 kpc. Such a large scale length was proposed by some authors 



( Ratnatu nga et al. 



2003 



Chang et al.l 



2011 



1989; Beers Sz Sommer-Larsen 1995; Chiba Sz Beer s 2000; Larsen Sz Humphreys 



rop 
3 



but it was always associated to a higher scale height, as showiM in Fig- 
ure El The observations indicate that the required thin and extended thick disk is very unlikely. In 
conclusion, varying the three parameters involved in the calculations within the ranges allowed by 
the literature does not solve the problem of the missing DM in the volume under analysis. 



4.2. Non-flat rotation curve 



The behavior of the Galactic rotation curve is still debated for R >11 kpc (e.g. jBinnev Sz Dehnen 



19971 ). but there is a general consensus that it is flat at the solar Galactocentric position. In- 
deed, recent observations find only a tiny negative gradient (IXue et al.l 120081 : iFuchs et al.l 120091 . 
dV /dR = —0.85 and —0.006 ± 0.016 km s _1 kpc -1 , respectively). Large deviations from a flat 
curve are therefore exc l uded, and theoretical mo de ls predict \dV JdR\ < 6 km s _1 kpc -1 (e.g., 



Dehnen Sz Binney 



1998 



Oiling Sz Merrifieldl l2000l ) . iLevine et al.1 (120061 ) also used a flat curve to 



measure the Galactic density distribution. 



5 The d a ta of 
OihaetalJ Jl999l) 



(2008) 



Juric et al 



Figure O are taken 


from Yamagata & Yoshii 


1992); Robin et al. 


Buser et al. (1998. 


19991); Siegel et al.l (2002): 


Cabrera-Lavers et al. 


(2008); de Jong et al.l J2OI0I); Carollo et al. (2010 


); Chane et al. (2011) 
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Fig. 2. — Absolute (left panels) and incremental (right panels) surface density calculated with differ- 
ent values of the thick disk scale length (upper panels) and scale height (lower panels), overplotted 
to the expectations of the models described in the text (grey curves) . 
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Fig. 3. — Measured thick disk scale length and height in the literature (full dots with lcr-error 
bars). The empty circle shows the value required to force the calculated to match the MIN 

model. 
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If the hypothesis (|VIII[) is dropped and a non-flat rotation curve is considered, the additional 



term 



1 



ttGR 



dV- 
8R 



Vdz 



(21) 



is added to the right hand term of both Equation [13] and [HI In the following, we will assume 
Vlsr=220 km s _1 , and the vertic al shear of the thi c k disk rotational velocity V(Z) will be repre- 
sented by the linear expression of lMoni Bidin et al.1 (|2012l ). that well approximates the underlying 
low-index power law in the Z-range under analysis. It is immediately evident from Equation [22 
that a decreasing rotation curve at R & (SX < km s _1 kpc -1 ) introduces a negative term to the 
calculation, thus increasing the discrepancy between the observations and the expectations of the 
DM halo models. On the contrary, a higher is obtained if the rotation velocity increases 

with R, and a certain amount of DM is thus allowed in the volume under analysis. However, 
^ = 10 km s _1 kpc -1 is required to match the minimum DM density deduced by the Galactic 
rotation curve (MIN model), and ^ = 16.5 km s _1 kpc -1 to match the SHM. Such steep rotation 
curves are excluded by observations. Moreover, we are comparing the observations with models 
whose density was fixed to return a flat rotation curve, hence we step into the contradiction that, 
while assuming a steeper curve to increase the measured mass, we simultaneously increase the 
amount of DM required by the models. Among the possibilities offered by the literature, the solu- 
tion closest to the expectations of a DM halo model is obtained assuming a ro tation curve inc r easing 



its steepness from 



z=o~- 



Oto 



i2L 

dR 



z=4kpc~ 7 km s _1 kpc -1 , as modeled by lKalberla et al.l (120071 ) 
for the Galactic gas. Even in this case, however, the MIN model is 2a higher than both the mea- 
sured £>i.5(Z) and S (4 kpc), while the SHM is excluded at the 5er level. In conclusion, assuming 
a non-flat rotation curve can alter the results, but insufficiently to justify the mismatch between 
the observations and the models. 



4.3. Alternative radial profiles of the dispersions 



A popular model for the radial dependence of cjtj, alternative to the hypothesis ([Xj) . is the 
assumption of a constant Toomr e Q-parameter. Its val i dity i s controversial, and it is disfavored 
by the numerical integrations of ICuddeford Amendtl (|1992l ). If the rotation curve is flat, the 
constant-Q model predicts that 

al{R) oc i? 2 exp(--^-) (22) 



h R, f 



(jAmendt Cuddefordl ll991). from which it immediately follows that 



<9o"u _ 2 

°" U "^ h R , p 



dR 



(23) 



When dropping the hypothesis ([Xj) to assume the constant-Q model, a choice about the radial 
dependence of Oy and UW must also be done. However, they have a lower incidence on the results, 



and from Equation! 



| it is easy to see that dafj /dR is almost indistinguishable in the two cases at 




12 3 4 

Z [kpc] 

Fig. 4. — Surface density as a function of Galactic height calculated assuming the radial constancy 
of the Toomre Q-parameter (full black curve) , and a linear radial decay of the dispersion (dashed 
curve). The dotted curves indicate the respective la strip. The models discussed in the text are 
also overplotted (grey lines). 
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the solar position, because R@ ~ 2h^ p . Hence, it makes little difference if they follow the trend 
of Equation [23l or a purely exponential decay. This is not the same for a 
derivative d 2 a 2 ] /dR 2 enters into Equation 1131 

Substituting the assumption (jX]) with Equation 
Equation HU with k\ replaced by 

15 6 6 



for (7jj 



u, because the second 



leads to a solution identical to 



R 2 



(24) 



The results are shown in Figured The curve of both increases its slope and shifts upward, and 
a certain amount of DM is allowed, but this is still insufficient when compared to the models: the 
expectation of the SHM is still ~ 4.5<r higher than both and £>i.5(Z), while the discrepancy 

reduces to 2a for the minimum-density MIN model. The assumption of the constant-Q model for 
the radial decay of crjj mitigates the gap between the observations and the models, but not enough 
to reconcile them. 



Neese &: Yosd (|1988l ) found that their measurements of cttj throughout the Galactic disk could 
be fitted with a linear radial decay, with slope —3.8 ± 0.6 km s^ 1 kpc -1 . The effects of drop- 



ping the assumption (|X|) in favor of linear radial profiles can be e xplored, altho u gh th ere is 
no theoretical support for this behavior. We will assume from Neese Sz Yoss ( 1988 ). and 
5^ = —6. km s -1 kpc~* from Casetti-Dinescu et al. ( 2011 ). whose results for a\j are identical 
to those of iNeese h Yossl ()1988l ). although they do not claim that the underlying trend should be 
linear. Fixing is more problematic due to the absence of measurements in the literature, but 
from the last two terms of Equation [13] it is easy to see that a more negative value decreases the 
slope of Ti(Z). Therefore, we will adopt the most favorable case [7W(-R)=constant. 

Assuming the radial linear trends given above, the calculation of from Equation 1131 is 

straightforward, and the results are shown in Figure [H The curve of has a much steeper 

slope than under the assumption (}Xj) . and S>i.5(Z) matches the expectation of the MIN model, 
agreeing with the SHM model within ~ la. However, still is offset by 10-15 M© pc -2 in the 

whole Z-range when compared to the MIN model, and 25-30 M© pc~ 2 with respect to the SHM 
curve. Thus, the theoretically unjustified linear radial profile of a\j still requires an additional 
ad-hoc correction, like a strong overestimate of the visible mass, or a thick disk scale height of 
hz,p=0-7 kpc, and it is therefore unlikely. 



4.4. A flared thick disk 



The flare of the Galactic thick disk is still an issue of de bate, and the observa tions have not 

(20 06|), for example 



provided unique evidence of its existence. ISiegel et al.l ((2002) and iDu et al 



argue in favor of a fl ared thick d i sk wit h sc ale height increasing w ith R, while ICabrera-Lavers et al 
(|2007l ). supported bv lBilir et al.l (|2008l ) and lYaz fc Karaalil (|2009i ). proposed a flare of opposite sign. 
Given these contradictory results, it is most probably safe to assume in our calculations that the 
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12 3 4 

Z [kpc] 

Fig. 5. — Surface density as a function of scale height when a flared disk is considered (black lines), 
for different values of the local flaring Kj/ = -jjf^- The models described in the text are also 
indicated with grey curves. 
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flare is negligible at the solar position, else a relevant variation of the thick disk scale height would 
have been detected beyond doubt. 



If we assume a flared thick disk, thus dropping the assumption pxp . the scale height in Equa- 
tion [9] becomes a function of R, and we have: 

SlnM_ i -I i + _4_.*5sy (25) 



9R ha, flare ^ ^R,p h Z,p( R ) dR 

The consequence of this new approach is the simple substitution of h^ p with hnji are in Equation[T3l 
plus the additional term 

j U OR ( hfiji are ) ^ ^ 

The new formulation, although analytically simple, hides many practical complexities: first, the 
expression for hz, P (R) must be accurate enough to provide a good approximation up to its second 
radial derivative. Additionally, the assumption (jX]) implicitly relies on the constancy of hz, P {R) 



(jvan der Kruit &: Searld Il98ll . [1982]) , and in a flared disk the radial behavior of the dispersions is 



harder to model. Here we will therefore limit to the first-order approxi mation that the flare is locall y 



small, hence well represented by a linear function of R (as proposed bv lCabrera-Lavers et al.ll2007l ). 
and not noticeably affecting the radial behavior of the dispersions. In this case, the additional term 
of Equation [26] is neglected, and a solution identical to Equation O is obtained, with k\ and k% 
respectively substituted by 

^ = 7r^ i^ + ^-(^--jr L -)' < 27 > 

-n-0 ' n RJlare <l R p tlR^p \ Kq liRJlare' 

k ' - 2 + 1 2 (28) 

, flare 

Rq 

The resulting profiles of the surface density are shown in Figure [5j The introduction of a flare in 
the formulation only changes the slope of the curve, that decreases in the case of a negative flare 
( d QR P < 0). A positive flare, however, does not return results that can be easily reconciled with the 
presence of a DM halo: £>i.5(Z) can be forced to match the MIN model assuming = 0.1, and 
the SHM model with d ^ R p = 0.16, but even in this case S(Z) still remains 15-20 M pc -2 lower 
than the expectations of this model in the whole z range. As in the case of a linear radial decay of 
the dispersion (Section 14. 3p . an additional correction to the formulation is therefore required. 



4.5. Alternative kinematical results 



We repeated the calculations replacing the kinematical measurements of iMoni Bidin et al 



(|2012i ) with other works in the literature, to check the incidence of the ass umed kinem a tics o n 



our results. Three previou s inv estigations a r e suit able for our purposes: iFuchs et al.l (|2009i ) 



Casetti-Dinescu et al.l (|201ll ). and lSmith et al.l (|2012l ). The results are shown in Figure [6l 
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Fig. 6. — Surfa ce density calculated b y use of the kinematical results of lFuchs et ajj ( 2009. short- 
dashe d curve), ICasetti-Dinescu et al.1 (|201ll . full curve with dotted la stripe), and ISmith et al. 
([2012], long-dashed curve). The expectations of the models discussed in the text are overplotted. 
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The data points of lFuchs et al.l (|2009l ) in the range Z=0-0.8 kpc were linearly fitted to derive 



the vertical profiles of cttj, cry, 0"w> an d UW to be inserted in the Equation (TIH). A quadrati c 
fit led to negligible differences. We assumed /iz iP =0.3 kpc and /ir iP =2.6 kpc (|Juric et al.l 120081 ). 
because their sample mostly c ompr ises thin disk stars. The curve of thus deriv ed is clearly 



unphy sical. However, iBinnevi (|2009l ) have already shown that the kinematical data of lFuchs et al 
(|2009l ) lead to inconsistent results when used to constrain the mass distribution in the Galaxy. 
The most likely cause of the problem is that their sample is a mixture of old thin and thick disk 
stars, with the incidence of the thick disk increasing with Z , caus ing too steep a vertical gradient 
of the dispersions. On the contrary, ICasetti-Dinescu et al.l (|201ll ) used a pure thick disk sample 
to measure the variation of the kinematics in the range Z=0.3-2.2 kpc. Unfortunately, they did 
not give the Z-profile of UW, but th eir measurement of the tilt angle a (strongly related to UW) 



agrees with iMoni Bidin et al.l (|2012l ). Hence we assumed the same profile of UW(Z) as before. 
In any case, the incidence o f this term is very limited at lower Galactic heights. The use of the 
Casetti-Dinescu et al.l (|201ll ) data returns the same general results previously discussed: well 
matches the expectations for the visible mass alone even in this case. Unfortunately, the results 
are weaker, both because of the larger observational errors, and because the expectations of the 
models are less distinct at lower Z. Thus, the less massive DM model is only la higher than the 
calculated solution. 



Smith et al.l (|2012l ) recently measured the change of the three-dimensional disk kinematics 
between Z=0.5 and 1.8 kpc, in three metallicity bins. The metal-rich group ([Fe/H]> —0.5) is 
probably a heterogeneous mix of thin disk sub-populations, and it is not suitable for our pur- 
poses. The metal-poor thick disk , dominating the sample with [Fe/H]< —0.8, has a unique spatial 
distribution (jCarollo et al.ll2010l ) poorly studied in the literature. Hence, we adopted the results 
for the intermediate-metallicity stars (—0.8 <[Fe/H|< —0.5), representative of a thick disk stellar 



population similar to that studied by lMoni Bidin et al.l (120121 ) . To derive the vertical trend of the 



kinematical quantities, we linearly fitted their results excluding their nearest bin (Z=0.7 kpc). We 
thus limited the incidence of a residual thin disk contamination, not efficiently excluded at low Z 
by the metallicity-only sample selection. Moreover, restricting the Z-range ensures that the linear 
fit is a good approximation of the underlying trend because close to the plane, as already discussed 
in Section [31 the vertical gradient of the dispersions is expected to vary rapidly with Z. We did 
not estimated an error on in this case, because the uncertainties associated to the kinemat- 



the solution obtainec 


with 


Smith et al. 


Casetti-Dinescu et al. 


f 


2011 


), again rou 



In conclusion, two more sets of independent kinematical results return results identical to ours 
when used in the calculation of S(Z), although with a lower significance. 
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Fig. 7. — Comparison between the surface density of the dynamical mass estimated by means of 
the Poisson equation (black curves), and of the mass enclosed between ±Z (grey curves), for the 
DM halo models SHM and N97. 
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4.6. DM halo detectability 

The failure to detect the presence of a DM halo naturally raises the question whether this 
feature is detectable by our method. In fact, the calculation does not measure the mass directly, but 



a vari ation of the gravitational potential, from which the mass density is derived. iMoni Bidin et al 



(|2010l ) noted that, if the DM halo was too uniform, extended, and poorly concentrated, it could 



cause a negligible change of the potential in the volume under analysis, thus resulting undetectable. 

The gravitational potential can be separated at any point as the sum of the potential of 
the dark and visible matter, <f> = 4>dm + '/'vis- It is easy to see that the dynamical surface density 
calculated through Equation (pQ) can therefore be expressed as the sum of two contribution, Ti(Z) = 
%4,dm(Z) + £<£,vis(-£), where 

^M(Z) = —[j o + (29) 

and analogously for y\s(Z). The potential of a spherical NFW DM distribution can be written 
as 

GMvMl + cs) 

0dm(s) = -g(c) , (30) 

r s c s 



where r s is a characteristic radial scale-length, c is the concentration par ameter (jLokas! l200lh 



s=r/(r s c) with r = V R 2 + Z 2 the radial spherical coordinate, My is the total mass within the 
Viral radius ry = r s c, and g(c) = (ln(l + c) — c/(l + c)) -1 . The two NFW models analyzed in 
this paper have M v = 5.2 x 10 11 M Q and c = 16.8 (SHM model), and M v = 6.85 x 10 11 M 
and c = 9 (N97 model). Inserting Equation ([30]) into Equation ([29]) . we can estimate X<^dm(Z), 
i.e. the dynamical DM mass that can be detected by means of the Poisson equation between ±Z 
from the Galactic plane. The results of this calculation for the SHM and N97 models are shown 
in Figure In both cases, the dynamical DM mass perfectly matches the DM mass expected by 
the models, indicating that all the mass of the DM halo enclosed between ±Z is detected by the 
use of Equation ([I]). It is therefore false that the small potential difference induced by the DM 
halo between and 4 kpc from the plane can induce an underestimate of the physical mass in the 
volume under analysis. 

It can also be argued that, even if EquationQ]should not fail, the change of the stellar kinematics 
between Z=l.h and 4 kpc induced by the presence of a DM halo could be too small to be detected 
by the kinematical measurements adopted in our work. To check this possibility, we compared the 
vertical gradient of the potential of a DM halo, fe^ , with the observed quantity 



dZ 



of*/ daf/ / 1 



YL-^-Tjw-[---^-). (31) 
obs hz,p oZ \R h^p) 



This is the vertical gradient of the gravitational potential estimated by means of the stellar kine- 
matics through Equation (|3|). The comparison, relative to the SHM model, is shown in Figure [HJ 
Clearly, the DM halo does not generate a strong vertical gradient of the potential: the surface den- 
sity of the expected DM mass enclosed within ±4 kpc is similar to the measured one (~60 M Q pc~ 2 ), 
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Fig. 8. — The vertical gradient of the potential caused by the SHM dark halo model (grey curve), 
and measured from stellar kinematics through Equation [31] (black solid curve), with its ±1<t strip 
(black dashed curves). 
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while is about a factor of two smaller than the observed potential gradient. Nevertheless, 

54 is measured with an error of the order of 10% (~200 km 2 s~ 2 kpc -1 ). This indicates that the 

dZ obs 

presence of a DM halo would had not passed undetected, because it would had caused a variation 
of the measured 54 about 4-5 times larger than the observational errors. In a forthcoming 

obs 

study, the orbit integration of a synthetic stellar population embedded in the Galactic potential 
will be used to test the reliablility of our method, and to analyze in more detail the causes of the 
systematics that can arise. 



4.7. Dark halo flattening 



So far, we have considered only spherical models for the DM halo, but the expected amount of 
DM in the volume under analysis changes noticeably if its shape is modified. The effects of varying 
the flattening of the models SHM and MIN are shown in Figure The models were calculated from 
Equations (fTUj) and (|2"0"j) . where R c, p c , and p p , dm were multiplied by 1/q to keep the resulting 
Galactic rotation curve the same (|011ingill995l ; lOlling fc Merrifieldll200ll ). Clearly, oblate models 
(q < 1) imply a larger quantity of DM in the solar vicinity, and they depart further from the 
observations. The expected surface density, on the other hand, decreases at increasing q, and 
highly prolate models approach the observational results. A perfect match requires a zero-density 
halo and is therefore never reached for any finite q, but a lower limit can be derived. We thus 
conclude that strongly prolate models are required, because within 2a (95% confidence level) we 
have q > 2 for the MIN model, and this lower limit is much higher for more common, higher-density 
models, for example 9 > 4 for the SHM. 



4.8. Other dark matter structures 



We have compared the observations with common models of DM halos, bu t other spatial 



distr ibution have been presented in the literature, like a DM disk and a ring (e.g., iKalberla et al 
20071 ). Nevertheless, these features were usually proposed in addition and not as an alternative to 
a spheroidal halo, thus enhancing the expected local DM density. 

The presence of a DM ring at R 3> Rq cannot be detected with our observations, but most of 
the mass of a DM disk would be included in the volume und er ana l ysis, a nd thus detectable by our 
method. The existence of such a feature, first proposed by iLakd (|1989l ). is of parti cular interest 



becau se it is a natural expectation of current ACDM models of Galactic formation (jPurcell et al 



2009). In 



igur e [TUI our results are compared to the expectations of t he massive DM disk model of 



Kalberlal ( 20031) . the thick, low-density model of I Pur cell et al.l (|2009l ) . and the thinner and denser 



ones from lRead et al.l (120081) . T he comparison shows that the dark disk proposed as an alternative 
to a DM halo (IKalberlal 120031 ) is excluded by the observations with very high significance (9er) 
and, as discussed in the previous sections, this mismatch cannot be corrected by simply altering 
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Fig. 9. — Comparison between the measured surface density (black line) and the model expectations 
(grey curves) with varying flattening parameter q. The dotted and dashed lines indicate the la 
and 3a stripes, respectively. Upper panels: MIN model. Lower panels: SHM model. 
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Fig. 10. — Comparison between the obse rved surface den sity (black line), and the e xpectations of 



DM d isk models from the literature. K07: iKalberlal (|2003l ) model; Pu-0.3 and Pu-0.15: iPurcell et al 



(|2009l ) models with local density pfT) i p/u=3 and 1.5 mM pc~ 3 , respectively, and scale height 
4.6 kpc; Re-1.0, Re-0.5, and Re-0.3: iRead et al.1 (|2008l ) models with p©,,dm=10j 5, and 3 mM pc -3 , 
respectively, and scale height 2.4 kpc. 
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one paramet er or one a s sumpt ion. On the contr a ry, th e local density of the DM disk models 
proposed by iRead et al.l (|2008l ) and iPurcell et al.l (|2009l ) is low, and they cannot be ruled out. 
Some constraints can be derived, though. The curves of both and E>i.s(Z) indicate that, 

within 2a, the thick DM disk of IPurcell et al.l (|2009l . scale height ^4. 6 kpc) mus t have a local 
density Pq,dm < 3 mM pc -3 , while the upper limit for the thinner IRead et al.l (|2008l ) model 
(scale height 2.4 kpc) is Pq,dm < 5 mM© pc -3 . However it must be taken into consideration that 
these results are sensitive to a change of a parameter or an assumption as discussed in Sections 14.11 
to 14.41 In any case, these DM disk models do not sustain the Galactic rotation curve, and their 
presence in addition to a DM halo would enhance the discrepancy between the observations and 
the expected DM density in the volume under study. 



5. DISCUSSION AND CONCLUSION 

The measurement of the mass surface density at the solar Galactocentric position between 1.5 
and 4 kpc from the Galactic plane accounts for the visible mass only. The DM density in the solar 
neighborhood, extrapolated from the observed curve of is Pq,dm = Oil mM^ pc~ 3 , at vari- 



ance with the general cons ensus that it must be in the range 5-13 mM pc (e.g. JWeber &: de Boer 



201Q : iGarbari et al.ll201ll ). Our recent measurements of the thick disk kinematics were used in the 
calculation, but the observed lack of DM is independent of this choice, because very similar results 
can be obtained by means of other kinematical results in the literature. The calculation relies on 
three input parameters and ten assumptions, but the observations cannot be reconciled with the 
DM halo modes modifying one of them. Altering them at will introduces enough freedom to force 
the solution to match the expectations of the most preferred model, but in this case an exotic 
series of unlikely hypotheses must be invoked. For example, a very thin thick disk (h^ jP =0.7 kpc), 
either very extended in the radial direction (h/j p =4.6 kpc) or strongly flared at the solar position 

Oh ' 

(-0^=0.1) make the solution coincide with the minimum DM local density deduced from the 
Galactic rotation curve, but the observational constraints exclude or disfavor these scenarios. On 
the other hand, the expected visible mass strikingly matches the observations without any effort, 
by use of the most probable assumptions. This coincidence lends weight to the interpretation of 
these results, because it is eas y to obtain an unphysic al solution if one or more wrong hypotheses 
are being made. For example, iMoni Bidin et al. showed that the results obtained under the 



assumption of a cross-term UW symmetric with respect to the Galactic plane, as an alternative to 
our hypothesis pvp . violate two minimum requirements: the surface density must at least account 
for the known visible matter, and it cannot decrease with Z. Thus, the excellent agreement be- 
tween the measured mass and the visible mass is unlikely to have been obtained by pure chance. 
We also demonstrated that our method cannot fail to detect the presence of a classical DM halo, 
because it causes a noticeable change in the stellar kinematics, one order of magnitude larger than 
the observational errors. 



The only viable solution to reconcile the observations with the models is the assumption of 
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a highly prolate DM halo, that can sustain a flat rotation curve with a negligible density at the 
solar position. Observational constraints on the DM halo shape are scarce a nd often controversial 



and, while spherical structures are usually preferred (e.g. Ilbata et al.l l200ll ; iMajewski et al 



Johnston et al 



2005 



Fellhauer et al.l 120061 ) . a prolate spheroid (q = 5/3) was invoked by 



2003 



Hclmi 



(|200J) and lLaw et al.l (|2005l ). Our observations suggest that q >2 is required even in the lower 
bound case of the least massive model, for agreement within 2c However, very prolate structures 
are atypical in cold dark matter simulations, which have problems in reproducing them (e.g., 
Dubinski Binnevl 119981 ), Thus, it must still be proven that a DM halo with a high flattening 



parameter is fully compatible with the current ACDM paradigm. 

A dark spheroidal component i s required to sustain t he Galactic rotation curve, observa tionally 



confirmed to be flat for R >5 kpc (|Kalberla et al 



2007 



Sofue et al. 



2008 



Xue et alJl200ah . In the 



presence of the visible mass only, Newtonian dynamics would predict a steep keplerian fall-off. Very 
noticeably, the calculation returns the unphysical result of a total surface density decreasing with 
Z if a rapidly decreasing curve is assumed. Hence, a flat rotation curve is required, while failing to 
detect the DM necessary to sustain it. This apparent contradiction is actually a confirmation that 
the calculation is reliable, because it is consistent with an observationally proven fact, although not 
with its expected explanation. 

In conclusion, the observations point to a lack of Galactic DM at the solar position, contrary 
to the expectations of all the current models of Galactic mass distribution. A DM distribution very 
different to what it is today accepted, such as a highly prolate DM halo, is required to reconcile 
the results with the DM paradigm. The interpretation of these results is thus not straightforward. 
We believe that they require further investigation and analysis, both on the observational and the 
theoretical side, to solve the problems they present. We feel that any attempt to further interpret 
and explain our results, beyond that presented in this paper, would be highly speculative at this 
stage. Future surveys, such as GAIA, will likely be crucial to move beyond this point. However, 
as our results currently stand, we stress that, while numerous experiments seek to directly detect 
the elusive DM particles, our observations suggest that their density may be negligible in the solar 
neighborhood. This conclusion does not depend on the cause of this lack of DM at the solar 
position. For example, if our results are interpreted as evidence of a highly prolate cold DM halo 
with q >2, this would have a local density lower than 2 mM pc -3 , i.e. more than a factor of four 
lower than what usually assumed in the interpretation of the results of these experiments. 

It is clear that the local surface density measured in our work, extrapolated to the rest of the 
Galaxy, cannot retain the Sun in a circular orbit at a speed of ~220 km s _1 . A deep missing mass 
problem is therefore evidenced by our observations. Indeed, we believe that our results do not solve 
any problem, but pose important, new ones. 
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